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Abstract 

Using the next-to-leading order QCD-corrected effective Hamiltonian, charmless exclusive non- 
leptonic decays of the B s meson into rj or rf are calculated within the generalized factorization 
approach. Nonfactorizable contributions, which can be parametrized in terms of the effective num- 
ber of colors iV| for PP and VP decay modes, are studied in two different schemes: (i) the one 
with the "homogeneous" structure in which (N° s )i (N% )2 (N^ s )iq is assumed, and (ii) 

the "heterogeneous" one in which the possibility of N^ S (V + A) ^ N^ S (V — A) is considered, where 
N% S (V + A) denotes the effective value of colors for the (V — A)(V + A) penguin operators and 
Nf{V - A) for the (V - A){V - A) ones. For processes depending on the iV^-stable cij such as 
B s — > (it, p)(jj' ,rj), the predicted branching ratios are not sensitive to the factorization approach we 
choose. While for the processes depending on the iV^-sensitive a« such as B s — > wr/ \ there is a 
wide range for the branching ratios depending on the choice of the N£ involved. We have included 
the QCD anomaly effect in our calculations and found that it is important for B s — > rf^rf\ The 
effect of the (cc) — > rj 1 mechanism is found to be tiny due to a possible CKM-suppression and the 
suppression in the decay constants except for the B s — > (f>rj decay within the "naive" factorization 
approach, where the internal W diagram is CKM-suppressed and the penguin contributions are 
compensated. 
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1. Motivation Stimulated by the recent observations of the large inclusive and exclu- 
sive rare B decays by the CLEO Collaboration [JXJ] , there are considerable interests in the 
charmless B meson decays @ . To explain the abnormally large branching ratio of the semi- 
inclusive process B — > rf + X, several mechanisms have been advocated p, |], [5], || and some 
tests of these mechanisms have been proposed 0. It is now generally believed that the QCD 
anomaly |3|, f|, [J plays a vital role. The understanding of the exclusive B — > rj'K, however, 
relies on several subtle points. First, the QCD anomaly does occur through the equation of 
motion || H when calculating the (S — P)(S + P) penguin operator and its effect is found 
to reduce the branching ratio. Second, the mechanism of cc — > rj', although proposed to be 
large and positive originally [fL0|, |TT| , is now preferred to be negative and smaller than before 



as implied by a recent theoretical recalculation [[12] and several phenomenological analyses 
P, Third, the running strange quark mass which appears in the calculation of the matrix 
elements of the (S — P)(S + P) penguin operator, the SU(3) breaking effect in the involved 
rj decay constants and the normalization of the B — > rf^ matrix element involved raise the 
branching ratio substantially. Finally, nonfactorizable contributions, which are parametrized 
by the N% s , gives the final answer for the largeness of exclusive B — > rj'K. ( We refer the 



reader to [14], |15[ for details 



It is very interesting to see the impacts of these subtleties mentioned above on the the 
exclusive charmless B s decays to an rf or rj. In addition to the essential and important QCD 



penguin contribution as discussed in |Tj| [17], it is found that the EW penguin contribution 
is important for some processes |18| , e.g. B s — > (n, p)(r), 4>) which the QCD penguin does 
not contribute to. However, the effects of QCD anamoly on the B s decay are not discussed 
in earlier papers. This motivates us to consider the contributions of anamoly effects in 
charmless B s decays. Another interesting topic we would like to study is the importance of 
the mechanism cc — > rj'. Besides, the running quark mass, the rf^ decay constant and the 
normalization of the matrix element involving rj^ are carefully taken care of in this Letter. 

2. Theoretical Framework We begin with a brief description of the theoretical framework. 
The relevant effective AB = 1 weak Hamiltonian is 

Q 10 

H eS (AB = 1) = -^\V ub V* (dO? + c 2 O u 2 ) + V ch V: q {c x O\ + c 2 2 c ) - V tb V*J2ciO<] + h.c, (1) 
V2 L i=3 

where q — d,s, and 

0\ = {ub) v _ A {qu) v _ A1 0\ = (cb) v _ A (qc) 

01 = (qb) v _ A (iiu) v _ A , O c 2 = (qb) v _ A (cc) 



V-A 5 
V-A 1 



v~a{v+a)i 
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°7(9) = ^(Qb) v _ A J2e q iqq') v+A ( v _ A ), O 8{10) = - {q a bp) v _ A ^ eg/^Jv+ACv-x). (2) 

q> q' 

with (q~iq2) v±A = ?i7^(l ±75)92- m Eq. (2), 3 _ 6 are QCD penguin operators and O7_i are 
electroweak penguin operators. Ci(fi) are the Wilson coefficients, which have been evaluated 



to the next-to-leading order (NLO) |T9|, |20| . One important feature of the NLO calcula- 



tion is the renormalization-scheme and -scale dependence of the Wilson coefficients (for a 
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review, see 



PHI). In order to ensure the [l and renormalization scheme independence for 
the physical amplitude, the matrix elements, which are evaluated under the factorization 
hypothesis, have to be computed in the same renormalization scheme and renormalized at 



the same scale as Ci(fi). However, as emphasized in [Ofl, the matrix element (0)f act is scale 
independent under the factorization approach and hence it cannot be identified with (O(fi)). 
Incorporating QCD and electroweak corrections to the four-quark operators, we can rede- 
fine Ci(fj,){Oi(n)) = c^ G (Oi)tree, so that df are renormaliztion scheme and scale independent. 
Then the factorization approximation is applied to the hadronic matrix elements of the op- 
erator O at the tree level. The numerical values for c!f ff are shown in the last column of Table 



I, where \x = mb(mb), A^- = 225 MeV, m t = 170 GeV and k 2 = ml/2 are used JEi 



MS 

In general, there are contributions from the nonfactorizable amplitudes. Because there is 
only one single form factor (or Lorentz scalar) involved in the decay amplitude of B (D) — > 
PP, PV decays (P: pseudoscalar meson, V: vector meson), the effects of nonfactorization 



can be lumped into the effective parameters a^ ff |22 



«* = + <&! + , a e 2 ?_i = d*_ x + c c i + X*-i) , (3) 

where c^ 2i _ 1 are the Wilson coefficients of the 4-quark operators, and nonfactorizable con- 
tributions are characterized by the parameters Xn and Xn-i- We can parametrize the non- 
factorizable contributions by defining an effective number of colors N* , called l/£ in 
as 1/Nf = (1/N C ) + x- Different factorization approach used in the literature can be 
classified by the effective number of colors Nf. The so-called "naive" factorization dis- 
cards all the nonfactorizable contributions and takes 1/Nf = 1/N C = 1/3, whereas the 



"large- N c improved" factorization \24\ drops out all the subleading 1/N C terms and takes 
l/Nf = 0. In principle, Nf can vary from channel to channel, as in the case of charm 
decay. However, in the energetic two-body B decays, Nf is expected to be process insen- 
sitive as supported by data p5] . If Nf is process independent, then we have a generalized 



factorization. In this paper, we will treat the nonfactorizable contributions with two dif- 
ferent phenomenological ways : (i) the one with "homogenous" structure, which assumes 
that [Nf)i ~ (^ ff )2 ~ • • • ~ (Nf) w j and (ii) the "heterogeneous" one, which consid- 
ers the possibility of Nf(V + A) ^ N* (V — A). The consideration of the "homogenous" 
nonfactorizable contributions, which is commonly used in the literature, has its advantage 
of simplicity. However, as argued in WQ, due to the different Dirac structure of the Fierz 



transformation, nonfactorizable effects in the matrix elements of (V — A)(V + A) operators 
are a priori different from that of (V — A)(V — A) operators, i.e. x^V + A) 7^ xty ~ A). 
Since 1/Nf = 1/N C + x > theoretically it is expected that 

Nf(V — A) = (Nf) i « (Nf ) 2 « (Nf) 3 « (Nf) 4 « (Nf) g « (Nf) w , 
Nf(V + A) = (Nf) s « (Nf) e Pa (Nf) ? « (Nf) s . (4) 

To illustrate the effect of the nonfactorizable contribution, we extrapolate N C (V — A) pa 2 
from B — > Dtt(p) [^] to charmless decays. The iV^-dependence of the effective parameters 
dj's are shown in Table I, from which we see that and a$ are iV^-stable, and the 

remaining ones are N* -sensitive. We would like to remark that while (13 and are both 
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Table I. Numerical values of effective coefficients a« at N% = 2, 3, 5,oo, where Nf = oo 
corresponds to af = cf . The entries for a3,...,aio have to be multiplied with fCP 4 . 





AT eff = 2 

C 


N cS = 3 

C 


N eS = 5 

c 


N eS = 00 




0.986 


1.04 


1.08 


1.15 


0,2 


0.25 


0.058 


-0.095 


-0.325 




-13.9 - 22.6z 


61 


120 + 18i 


211 + 45.3i 


CI4 


-344- 113i 


-380 - 120i 


-410 - 127i 


-450 - 136i 


a 5 


-146 - 22.6i 


-52.7 


22 + 18i 


134 + 45.3i 


a 6 


-493 - 113z 


-515 - 121i 


-530 - 127i 


-560 - 136i 


07 


0.04 - 2.73i 


-0.7-2.73z 


-1.24- 2.73i 


-2.04 - 2.73i 


a 8 


2.98- 1.37? 


3.32 -0.9i 


3.59 - 0.55z 
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a 9 


-87.9 - 2.73i 


-91.1 -2.73i 


-93.7- 2.73i 


-97.6 - 2.73i 


a 10 


-29.3 - 1.37* 


-13 - 0.91i 


-0.04- 0.55?; 


19.48 



Nf- sensitive, the combination of (03 — a 5 ) is rather stable under the variation of the Nf 
within the "homogeneous" picture and is still sensitive to the factorization approach taken 
in the "heterogeneous" scheme. This is the main difference between the "homogeneous" and 
"heterogeneous" approaches. While a 7 , a 8 can be neglected, a 3 , a 5 and a w have some effects 
on the relevant processes depending on the choice of Nf. 

Before carrying out the phenomenological analysis, we would like to discuss the dynamical 
mechanism involved. We first come to the QCD anomaly effect. As pointed out in H, |HS| ], 
the QCD anamoly appears through the equation of motion 

d^sjtfss) = 2m s si l5 s + -^G^GT. (5) 

47T 

Neglecting the u and d quark masses in the equations of motion leads to 

a* 



2 

mi 



(rf\-£GG\0) = ftm*, (6) 

and hence 

(v'Mo)^-i^- (7) 

To determine the decay constant fi, we need to know the wave functions of the physical 7/ 
and rj states which are related to that of the SU(3) singlet state rj and octet state i] 8 by 

r]' = r/s sin 6 + 770 cos 9, 77 = r/s cos 9 — r? sin 9, (8) 

with 1 770) = ~^\uu + dd + ss), \r] 8 ) = -^\uu + dd — 2ss) and 9 ~ —20°. When the r] — r/ 
mixing angle is —19.5°, the 7/ and 77 wave functions have simple expressions [FJ: 

\rf) = —j=\uu + dd + 2ss), \rj) = —j=\uu + dd — ss). (9) 
V6 ' ' a/3 
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At this specific mixing angle, ft = ^ft, in the SU(3) limit. Introducing the decay constants 
/ 8 and f by 

(0\Al\ Vo ) = if oPfi , (0\Al\ Vs ) = if sPll (10) 
then ft and /*/ are related to / 8 and f by |31 



ft = JjL sin6 + JjL cos 0} ft = _ 2 -| S in 9 + A cos 9. (11) 

Likewise, for the rj meson 

ft = _| cos # _ _| sin Bj f s = _ 2 _| cos ^ - A sin 6 . (12) 

From a recent analysis of the data of r], rf — > 77 and 77, rj' — > 7177 [j32]l , / 8 ( ) and 9 have been 
determined to be 

j- = 1.38 ±0.22, A = 1.06 ± 0.03, 9 = -(22.0 ± 3.3)°, (13) 

In In 

which lead to 

ft = 99MeV, ft = -108MeV, ft = 47MeV, ft = 131MeV. (14) 



For the u and d quarks involved, we follow mM to use 



(r]'\wy 5 u\0) = (r]'\dj 5 d\0) = ry (r]'\sj 5 s\0), (15) 
with r (i) being given by 





/2/o 2 


-/I 


cos 6> + sin 6* 




JVi 


-/o 2 


cos — v2 sin 6* ' 


1 


\/2/o 2 


-fi 


cos 6 1 — v^2 sin 9 


2 


\/2/8 2 


-fi 


cos + i sin 9 



y «/ U ./ 8 WO 1/ — V ^ 0111 u flf!\ 

We next discuss the cc — > r/^ mechanism. This new internal J^-emission contribution 
will be important when the mixing angle involved is V c bV* s , which is as large as that of 
the penguin amplitude and yet its effective parameter a| is larger than that of penguin 
operators. The decay constant /i/ , defined as (0| 07^750^') = ifjj° q^, has been determined 
from the theoretical calculations [[II], [H], [TIJ and from the phenomenological analysis of the 
data of J/ij} — > rjc'-f, J/ip — > 7/7 and of the 777 and 7/7 transition form factors || [13| . In 
the presence of the charm content in the 770, an additional mixing angle 9 C is needed to be 
introduced: 

|?7o) = —j= cos9 c \uu + dd + ss) + sm9 c \cc), 



V3 

1 

7s 



\r] c ) = -= sin 9 C \ uu + dd + ss) + cos 9 C \ cc). (17) 
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Then f5 = cos 6 tan 6 c f Vc and f° = — sin ^ tan ^c/^, where the decay constant f Vc can be 
extracted from r\ c — > 77, and # c from J/-0 — > rj c ^ and J/ip 7/7 [f|. In the present paper 
we shall use 

= -6 MeV, ^ = - tan 6f', = -2.4 MeV, (18) 
for 9 = —22°, which are very close to the values 

/c = -(6.3 ± 0.6) MeV, /° = -(2.4 ± 0.2) MeV (19) 



obtained in [13 



In the following we will show the input parameters we used. One of the important 
parameters is the running quark mass which appears in the matrix elements of (S — P)(S+P) 
penguin operators through the use of equations of motion. The running quark mass should 
be applied at the scale p ~ mj because the energy release in the energetic two-body charmless 



decays of the B meson is of order m b . In this paper, we use [£7 



m u (m b ) = 3.2 MeV, m d (m b ) = 6.4 MeV, m s (m b ) = 105 MeV, 

m c (m b ) = 0.95 GeV, m b (m b ) = 4.34 GeV, (20) 

in ensuing calculation, where we have applied m s = 150 MeV at // = 1 GeV. 

It is convenient to parametrize the quark mixing matrix in terms of the Wolfenstein 



parameters: A, A, p and i] , where A = 0.81 and A = 0.22 [28]. A recent analysis of all 



available experimental constraints imposed on the Wolfenstein parameters yields |29| 

p= 0.156 ±0.090, fj = 0.328 ± 0.054, (21) 

where p = p(l — 4r) and f\ = r](l — 4r), and it implies that the negative p region is excluded 
at 93% C.L.. In this paper, we employ the representative values: p = 0.16 and r\ = 0.34, 
which satisfies the constraint \/ p 1 + rj 1 = 0.37. 

Under the factorization approach, the decay amplitudes are expressed as the products of 
the decay constants and the form factors. We use the standard parametrization for decay 
constants and form factors | 23j| . For values of the decay constants, we use f n = 132 MeV, 



f K = 160 MeV, f p = 210 MeV, f K * = 221 MeV, f u = 195 MeV and = 237 MeV. 
Concerning the heavy-to-light mesonic form factors, we will use the results evaluated in the 
relativistic quark model |23|, |3(| by directly calculating B s i u ,d) P and B s ( u ,d) ~^ V form 



factors at time-like momentum transfer. Denoting r/ s = ss, the explicit values for the form 
factors involved are F^ Vs {0) = 0.48, fJ'^(O) = 0.44, and A% aK * (0) = 0.28, which are larger 
than BSW model's results fTE| . The q 2 dependence of the matrix element, parametrized 
under the pole dominance ansatz, are found to have a dipole behaviour for A ,Fi, and a 
monopole one for F . In the following, we will use the exact value calculated at the relevant 
kinematical point in this paper. Note that these matrix elements should be used with a 
correct normalization fl4j], for which to a good approximation, we take Fi q V (0) = : ^=F 1 q Vs (0) 

andF^'(0) = ^F 1 B ^(0). 
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3. Phenomenology We are now ready to discuss the phenomenology of exclusive charm- 
less rare B s decays. To illustrate the issue of N£ s dependence ( which means the different 
factorization approach) of theoretical predictions, we will begin with B s — > (7r, p, u)rf\ Un- 
like the case of B — > (n,p)r]('\ B s — > (n,p,u})r)^ do not receive the anomaly contribution 
from the (S — P)(S + P) penquin operators due to the particle content of B s and n(p, uj). 
The decay amplitude for B s ° — > nrj^ reads 



A(B S - v (\) 



V ub V*a 2 - V tb V t 



ts 



-(-07 + 09) 



>X (B s /),.) ; 



(22) 



where 



X^ (,) -) = (A(uu) v _ A \0)( V ^\(sb\ 



JB S ) = -i^{ml s -ml (l) )F^\rnl). (23) 



Since the internal W-emission is CKM-suppressed and the QCD penguins are canceled out 
in these decay modes, B s — > ir(p)r]^ are dominated by the EW penguin diagram. The 
dominant EW penguin contribution proportional to ag is iV£ ff -stable, whereas the internal 
W contribution a 2 is - sensitive. Within the "heterogeneous" nonfactorizable picture, a 2 
is fixed and thus the predicted branching ratio is rather stable under the variation of N* s 
as shown in the last four columns in Table II. However, a 2 varies within the "homogeneous" 
nonfactorizable scheme and thus the predicted branching ratios do show a Nf* dependence. 
We would like to emphasize that although B s — > (it, p)i]^ are dominated by the EW penguin 
diagram, the internal W diagram does make some contributions to this decay mode. Since a 2 
changes sign from iV? ff = 2, 3 to N° s = 5, 00, the interference pattern between the internal W 
diagram and the EW penguin diagram will change from the destructive to the constructive 
one. It is thus easy to see that for Nf 1 = 2 there is a larger destructive interference between 
the internal W diagram and the EW penguin contribution and the predicted branching 
ratio is the smallest one among the first four columns in Table II, whereas for N° s = 00 
constructive interference takes the role and the branching ratio increases. 

While QCD penguin diagrams are canceled out in B s — > (ir,p)r]('\ B s — > uirj^ gets 
enhanced from the QCD penguin diagram. The decay amplitude for B s — > corj^ is 



Gf 
V2 



V ub V*a 2 - V tb V t 



ts 



2(a 3 + o 5 ) + -(a 7 + a 9 ) 



(24) 



where 



* (iWV) = ( U \(vu) v _ A \0)( V O\^) v _ A \B a ) 



V^f^F^imDie-psJ. (25) 



From Table II, we see that there is a wide range of predictions for the branching ratios. 
This process is QCD penguin dominated, except when the "naive" factorization is used or 
(Nf(V - A),Nf(V + A)) = (2,5) where there are large cancellations between the QCD 
penguin contributions ( i.e. a 3 + a 5 ). The largest branching ration predicted for B s — > urf ) 
occurs when we use the "large- N c improved" factorization, where the EW penguin and QCD 
penguin have constructive interference. 
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Next, we discuss B s — > rj'K decay, which has the decay amplitude 



]-a w + (2a 6 - a 8 



v tb v td 



rn 



04 



K 



(m s + m d )(m b 



+ (a 3 - a 5 - a 7 + a 9 ) Xf* 1 ^ > + (a 3 - a 5 + -a 7 
+ (a 3 -a 5 -a 7 + a 9 )Xf^') 
a 3 + a 4 




) 



+ 



+ Oe - -as] 



111 

05 H CI7 — — Og — — Clio 




m s (m b - m s ) 



where 



(K°\(sd\ 

{v\(qq) v - 



(B s K,r) ) 



(26) 



JO)(r ] '\(db) v _ A \B° s ) = - i f K (m 2 Bs -m 
\0}(K\sb) v _jB° s ) = -zfUml s - m[ 



-lB 3 r) 



-B S K 



2 \ 

m 



n 



(27) 

Due to the QCD anamoly, there is an extra (ft / ft — 1) term multiplied with a 6 in X d BsK,ri \ 
whose presence is necessary in order to be consistent with the chiral-limit behaviour of the 
(S-P)(S+P) penguin matrix elements \T%\. Though penguin diagrams play the dominant 
role, the internal W diagram and the mechanism of the cc pair into the rf do have some 
nonneglible effects when N* s = 2 and = 00 where a 2 gets the larger values. Due to the 
large cancellation, the EW penguin has only tiny effect and can be neglected. The monotonic 
decrease of the branching ratio from N° s = 00 to N% = 2 within the "homogeneous" 
nonfactorizable picture can be understood from the behaviour of the QCD penguin i.e. 
the destructive interference between 0(3,5) an d «(4,6) : as N* s decreases, 0(3,5) contributions 
increase and hence the branching ratios decrease. As we already mentioned before, a 3 and 
05 are N% -sensitive while a 3 — as is stable under the variation of N^r and then the predicted 
branching ratio is iV^-stable within the "homogeneous" factorization approach. 

There exist some general rules for the derivation of the formula from B — > P a P b to 
its corresponding B — > V a P b and B — > P a V b . These general rules can be written as: (i) 
For X( BPa ' Pb ^ to X^ BVa,Pb \ replace the term wip /[(mi + m 2 ){m 3 — m.4)] by — mp /[(mi + 
m 2 )(m 3 +m4j] and the index P a by V a , (ii) discard the (S — P)(S + P) contribution associated 
with x( BsPa ' Vb * > and 0(5,7) ~~ > — a (5,7) if they contribute. Thus, the factorizable amplitude of 
B s — > rj'K* can be readily obtained from the B s — > rj'K one and reads 



A(B S ^r,K 



v tb v; d 



V u bV* d a 2 X^ 



(a 4 - -a w )X 



-(B a K*, v ) 



V cb V: d a 2 X^ 



(BsV ,K*) 



1 m 1 

(a 6 - ^a 8 ) — - — \ 

2 m s [m b + m 




a 4 - a 5 + -07 - -a 9 



(B s K*,r, ) 



rOio 
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+ (a s -a 5 -a 7 + a 9 ) X^ K "^ + (a 3 - a 5 + ^a 7 - Ja 9 )Xf 



+ («3 



a 5 - a 7 + a 9 )Xf^ (,) ) |, 



(28) 



with 



1 (su) v _ A \0) (/> | («6) v _ A \B S ) = 2f K ,m K ,F^> {,) (m 2 K ,)(e -p Bs ), 



( V (\qq) v _ A \0)(K*-\(sb) v _ A \Bj 



BsK* , 



It is interesting to see that since there is no term in penguin contribution is 

reduced substantially and so does the branching ratio. With a reduced penguin contribution, 
the involved internal W diagram and the mechanism of the cc pair into the 7/ become more 
important than those of the B s — > Kr( . The larger branching ratios in columns denoted by 
I a and II b are due to the constructive interference between the internal W diagram and the 
penguin contribution. 



We are now coming to the most complicated process B 
amplitude 



A(B S 



rjrj 



y rjrj' , which has the decay 
W£ (a2*J W) + a 2 Xf ^) + V cb V c * s (oa^W) + a 2 Xf^)) 



- v tb v* 



1 1 1 

a 3 + a 4 — a 5 + -a 7 — -a 9 — -aio 



1 m| 

2 m s [m b — m s 



, f S J) S 



+ (203 - 2a 5 - ia 7 + ^a 9 ) Xf*™') + (a 3 - a 5 - a 7 + a 9 )Xf 



+ 
+ K 



111 

a 3 + a 4 - a 5 + -a 7 - -a 9 - -ai 



2 m s (m b — m s 



f r) , 



\ X {BW,r,) 



+ (2a 3 - 2a 5 - ^a 7 + ^a 9 ) Xf^'^ + (a 3 - a 5 - a 7 + a 9 )Xf 



,»?) 



(30) 



with 



= W\m v -M(v\(sb) v _ A \B s ) = -ifi{m%.-ml)F*«(rnZ,), 
X?'™ = (v\m v . A \0)W\(sb) v . A \B s ) = -if^ml s -ml)F B ^(m^). 

The destructive interference between X^ffJ'V and X^f^^ makes the internal W^-emission , 
cc — > 77 ^ and the corresponding penguin contributions smaller. The EW penguin is smaller 
than the QCD penguin by an order of magnitude at the amplitude level and hence can be ne- 
glected. The dominant QCD penguin contributions are governed by Xf*™') an d xf^'^\ 
which have the constructive interference. The (a 3 — as) term is positive, contrary to the 
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negative 04 and terms, and becomes smaller when N° s increases within the "homoge- 
neous" nonfactorizable structure. Thus a monotonic increase of the branching ratio when 
Nf* increases comes mainly from this reduced destructive interference, within the "homoge- 
neous" nonfactorizable structure. Similar arguments are also applied to the "heterogeneous" 
structure. 

With the general rules (i) and (ii) mentioned before, the decay amplitude for B s — > <jnf) 
can be easily obtained from B s — > 7777': 



' 111 

03 + 04 — a$ + -07 — —ag — 2 Ol 



- v tb v* 



° 6 m s (m b + m s ) 

+ (2a 3 - 2a 5 - ^a 7 + ^a 9 ) J^*^ + (a 3 - a 5 - 

+ ( a 3 + a 4 + a 5 - ^a 7 - ^9 - ^10 





with 



While the internal W diagram is subject to the CKM-suppression, the cc — > rf^ mechanism 
suffers from the suppression in the decay constant and thus B s — > (jnf ) is dominated by the 
penguin contribution. Due to the cancellation among the different a^i = 3,4,5,6)'s, effect 
of the QCD penguin, though still dominant, are reduced substantially. Within the "homo- 
geneous" nonfactorizable picture, we find a monotonic decrease of the branching ratios when 
N% increases, which comes from a monotonic decrease of the QCD penguin contributions as 
N* s increases. Since the QCD penguin contributions are reduced, the EW penguin contribu- 
tions become important. It is found that the interference pattern between the QCD and EW 
penguin is destructive except for the "large- N c improved" factorization approach, where a 
constructive interference exists and a very dramatically suppressed QCD penguin contribu- 
tion appears. The strength of the destructive interference depends on N% and its effect is to 
reduce the QCD penguin contribution without changing the trend of reduced penguin when 
Nf 1 increases. A dramatic destructive interference among the penguin contributions occurs 
for B s -> 0/) when Nf{V-A) = Nf{V+A) = 00 and {Nf{V-A), Nf{V+A)) = (2, 00), 
thus the internal W diagram and the cc — > rf^ mechanism contribution becomes important 
relative to the other cases and the branching ratio is the smallest in these situations. 



4. Summary and Discussions To summarize, we have studied charmless exclusive non- 
leptonic B s meson decay into an 77 or r( within the generalized factorization approach. Non- 
factorizable contributions are parametrized in terms of the effective number of colors N* s 
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Table II. Average branching ratios (in units of 10 6 ) for charmless B s decays to rf and r\. Predictions 
are for k 2 = ml/2, r\ = 0.34, p = 0.16. I denotes the "homogeneous" nonfactorizable contributions 
i.e. Nf(V - A) = Nf(V + A) and (a,b,c,d) represent the cases for A c cff =(oo,5,3,2). II denotes 
the "heterogeneous" nonfactorizable contributions, i.e. N^ S (V — A) / N^ S (V + A) and (a' ,b' ,c') 
represent the cases for N^ S (V + A) = (3,5,oo), where we have fixed N* S (V — A)=2 (see the 
text) 



Decay 


/a 


h 


h 


h 




lib' 


Hcf 


B s -» 7TT?' 


0.25 


0.17 


0.13 


0.11 


0.11 


0.11 


0.10 


B s — > irq 


0.16 


0.11 


0.08 


0.07 


0.07 


0.068 


0.067 


B s - PV ' 


0.70 


0.47 


0.36 


0.30 


0.30 


0.30 


0.31 


B s — > pi] 


0.45 


0.30 


0.24 


0.19 


0.19 


0.19 


0.20 


B s — > LOT]' 


6.9 


0.9 


0.012 


2.14 


0.48 


0.03 


0.83 


B s — > a;?] 


4.45 


0.63 


0.008 


1.39 


0.31 


0.02 


0.54 


B s - r/X° 


1.25 


1.07 


1.01 


1.00 


1.27 


1.51 


1.90 


B s -> r/if° 


1.35 


0.81 


0.68 


0.76 


0.75 


0.74 


0.72 


B s -» r/ii *° 


0.49 


0.35 


0.32 


0.26 


0.49 


0.60 


0.80 


B s - r/X*° 


0.45 


0.05 


0.02 


0.24 


0.24 


0.24 


0.25 




47.4 


41.8 


38.3 


34.4 


39.5 


44.1 


51.5 


s s - 


26.6 


24.9 


23.8 


22.4 


33.8 


43.9 


62.2 


B s -> 7777 


20.3 


17.1 


15.1 


12.8 


11.6 


10.7 


9.1 




0.44 


0.59 


2.29 


6.20 


4.41 


3.11 


1.66 


B s -> 


0.04 


0.91 


2.29 


4.92 


2.28 


0.92 


0.10 



and predictions using different factorization approaches are shown with the N* dependence. 
It is found that for processes depending on the iV? ff -stable a,'s such as B s — > (ir,p)r]('\ the 
branching ratios are not sensitive to the factorization approach we used. While for the pro- 
cesses depending on the iV^-sensitive a^'s such as the B s — > ur}('\ the predicted branching 
ratios have a wide range depending on the choice of the factorization approach. The effect of 
the QCD anomaly, which is not discussed in the earlier literature, is found to be important 
for the B s — > rj^rj^. We also found that the mechanism (cc) — > rj' , in general, has smaller 
effects due to a possible CKM-suppression and the suppression in the decay constants except 
for the B s — > tpi] under the "large- N c improved" factorization approach, where the internal 
W diagram is CKM-suppressed and the penguin contributions are compensated. 

In this Letter, we, following the standard approach, have neglected the W-exchange and 
the space-like penguin contributions. Another major source of uncertainties comes from the 
form factors we used, which are larger than the BSW model's calculations. Although the 
Wolfenstein parameter p ranges from the negative region to the positive one, we have "fixed" 
it to some representative values. The interference pattern between the internal W diagram 
and the penguin contributions will change when we take a different sign of p. We will study 
these form factor- and CKM- dependence involved and all the B s — > PP, VP, VV decay 
modes in a separate publication. 
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